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Abstract. The starburst galaxy M82 has been mapped 
at A 1.2 mm using the IRAM 30-m telescope and the 
MPIfR 19-channel bolometer. The influence of the side- 
lobe pattern and the 12 CO(2— 1) line has been carefully 
analyzed. Based on this analysis, we conclude that the 
cold dust does not extend far into the halo as has been 
previously claimed. The total mass of dust in the inner 
3 kpc of the galaxy is 7.5 • 10 6 M Q . Assuming a solar gas- 
to-dust ratio, the inferred total mass of gas is 7.5 • 10 8 M Q 

Unlike in any other galaxy we find the CO(2— 1) line 
emission to be clearly more extended than the A1.2mm 
dust emission. A likely explanation for this exceptional 
behaviour, seen for the first time in galaxies, is heating of 
molecular gas by low-energy cosmic rays and soft X-rays 
in the halo of M82. 

Key words: dust - galaxies: individual: M82 - galaxies: 
ISM - galaxies: starburst - radio continuum: galaxies - 
radio continuum: ISM 



1. Introduction 

The prototypical starburst galaxy M 82 is the prime exam- 
ple of a galaxy in which the violent star formation activity 
gives rise to the formation of a bipolar outflow and an asso- 
ciated extended halo, which is visible in various regimes of 
the electro-magnetic spectrum. Outflowing material was 
first evident in the Ha light (Lynds & Sandage 1963, 
McCarthy et al. 1987, Bland & Tully 1988). The kine- 
matics of the outflowing material and the geometry of the 
cone have been thoroughly worked out by McKeith et al. 
(1995). The recent detection of Ha emission ~11 kpc away 
from the plane of M 82 (Devine & Bally 1999) determines 
the cap of this outflow. The whole scenario is further cor- 
roborated by the existence of vertical magnetic field lines 
(Reuter et al. 1994), along which relativistic particles par- 
take in the outflow, thus forming a radio halo (Seaquist & 
Odegard 1991). This halo exhibits a filamentary structure 
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away from the plane (Reuter et al. 1992). Hot gas has been 
also found in some kind of halo which extends for several 
arcminutes along the minor axis (Fabbiano 1988, Schaaf 
et al. 1989, Bregman et al. 1995). This X-ray emission cor- 
relates well with the Ha if observed with sufficiently high 
angular resolution (Watson etal. 1984). 

In this context, the search for neutral gas and dust 
away from the plane of M82 has seen a number of at- 
tempts. Owing to the close interaction between M81 
and M82, neutral hydrogen is seen enveloping M82, but 
whether this gas stems from the M 81 or has been expelled 
from M82 is still a matter of debate (Yun etal. 1994). 
Molecular gas associated with the dusty filaments outside 
the plane of M 82 has been reported by Stark & Carlson 
(1984) and has been observed by Sofue etal. (1992) out 
to a projected distance of ±2 kpc. Clear evidence for the 
existence of scattering dust away from the plane of M 82 
comes from measurements of optical polarization (Bing- 
ham etal. 1976, Notni etal. 1981, Neininger etal. 1990). 

Owing to its high luminosity in all spectral bands, 
M 82 has been a prime target for first-light experiments 
in the mm- and submm-regime (Elias et al. 1978, Jura 
etal. 1978). With the improvement of bolometric mea- 
surements detailed studies of the distribution of the cold 
dust in this galaxy have become feasible; thus the recent 
past has seen an increasing number of such investigations 
(A2 mm: Kuno & Matsuo 1997; Al mm: Thronson et al. 
1989, Hughes etal. 1990, Kriigel etal. 1990a, Alton etal. 
1999; AO. 5 mm: Jaffe etal. 1984, Smith etal. 1990, Alton 
etal. 1999). Carlstrom & Kronberg (1991) have shown 
that at A3 mm thermal free-free emission still dominates 
the overall spectrum. 

Apart from these dedicated bolometric measurements, 
continuum maps at millimeter wavelengths have been pro- 
duced as by-products in various observations of spectral 
lines (e.g. Neininger etal. 1998). As a result of the discov- 
eries of a galactic wind and halo in M82, the more recent 
bolometric measurements aimed at detecting a dust halo 
at mm wavelengths. CO emission away from the plane 
was reported by Stark & Carlson (1984), Nakai etal. 
(1986) and by Sofue et al. (1992). An outflow of molecular 
gas was claimed by Nakai etal. (1987). A dust continuum 
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halo was first mentioned by Hughes etal. (1990). Obser- 
vations with higher resolution (Kuno & Matsuo 1997) 
suggest emission at A2 mm out to 400 pc from the plane. 
Recently Alton etal. (1999) reported a dust outflow from 
the central region of M82, based on their sub-mm images. 

Here we report observations of M82 in the A1.2 mm 
continuum performed with the MPIfR 19-channel bolome- 
ter at the IRAM 30-m telescope. In Sect. ||we describe the 
observations and data analysis, with consideration of the 
possible influence of the error beam. In Sect. || the distri- 
bution of the cold dust in M 82 will be presented, with a 
discussion conducted in Sect. |], along with a comparison 
with other published measurements. In Sect. |^ we give a 
short summary of our results. 

2. Observations and data reduction 

2.1. Observations 

We used the 19-channel bolometer of the MPIfR (Kreysa 
et al. 1993) in the Nasmyth focus of the IRAM 30-m tele- 
scope during the spring session in 1997. The weather con- 
ditions were fair during our night-time runs, with rather 
stable zenith opacities between 0.2 and 0.3 at 230 GHz. 

The 19 channels of the bolometer are arranged in a 
closely packed hexagonal array, with beamsizes of 11" 
FWHM and spacings of 20". The calibration was per- 
formed by mapping Uranus every morning. The point- 
ing and focus checks were made at regular intervals dur- 
ing the observations using nearby quasars (in particular 
1308+326). The map was centered on the near-infrared 
peak at a 50 = 09 h 51 m 43'/4, A 50 = 69°55'01" (Joy etal. 
1987). The data were taken in the standard mapping mode 
where the field is scanned at a constant speed of 4" per 
second in azimuth, with subsequent scans being spaced 
by 4" in elevation. The map size was chosen large enough 
(320" x 270") to assure a good baseline determination in 
the presence of extended emission. We thus made sure that 
the field was mapped sufficiently far out perpendicular to 
the major axis of M82, such as to trace the full extent of 
any dust component in the halo. While scanning the sub- 
reflector was wobbled at a 2 Hz rate with an amplitude of 
45", which is small enough to ensure a stable beam pat- 
tern. Due to the small wobbler throw the OFF position 
was not necessarily free of emission, but using the algo- 
rithm of Emerson et al. (1979) this should have no effect 
on the resulting map. We also mapped Mars once and the 
point source 3C 273 twice each in order to obtain a de- 
tailed beam pattern. In addition, more maps of Mars and 
other calibration sources were used to properly check the 
calibration. 

2.2. Standard data reduction 

The whole field was covered six times. For each coverage 
we performed the data reduction separately using NIC 



(GILDAS software package), starting with the subtraction 
of a third-order baseline. 

The scanning mode described in the previous section 
leads to maps which contain a superposition of a positive 
and a negative image of the source (double-beam maps). 
In principle the deconvolution can be done by dividing the 
Fourier transform (FT) of the measured intensity distribu- 
tion by the FT of the wobble function and transforming 
back to image space. As the FT of the wobble function 
is a sine wave with zeroes at the origin and at harmon- 
ics of the inverse wobbler throw, it would cause problems 
at these spatial frequencies. To avoid these problems we 
used the algorithm of Emerson et al. (1979) to restore the 
double-beam maps into equivalent single-beam maps. 

Thereafter we combined the 19 channels and calibrated 
the resulting map with respect to Mars and Uranus (as- 
suming brightness temperatures of 198.5 K and 97.2 K 
respectively). In a last step the four best coverages were 
averaged and regridded onto equatorial (B1950) coordi- 
nates. (Two coverages were ignored because of bad S/N 
ratios.) This final map has a sensitivity of 5 mJy per beam 
and a angular resolution of 12" (HPBW). 

2.3. Influence of the error beam 

In order to check the influence of the error beam on the 
extended emission we applied a CLEAN algorithm de- 
veloped for single-dish maps (Klein & Mack 1995) to 
each coverage using a radially symmetric antenna pattern. 
This pattern was obtained in the following way: A field of 
4'5 x 3'5 size was mapped centred on 3C273. This an- 
tenna pattern was then rotated in small steps (5°), and 
these individual (rotated) maps were averaged to yield a 
nearly radially symmetric pattern. The subtraction of the 
CLEANed map from the original one showed that the er- 
ror beam contributes less than 25% of the rms noise level 
to the extended emission (> 20" away from the central 
position). 

Even the sidelobes can be ignored in our case, because 
at A 1.2 mm they are of the order of 2%, corresponding 
to roughly 5mJy/beam in the M82 maps. Averaging four 
coverages with different parallactic angles (and accord- 
ingly different positions of the sidelobes) the effect of the 
antenna pattern on the extended emission falls well below 
the rms noise of our final map. 

This extra analysis make sure that the extended emis- 
sion seen in our map is real and not affected by the error 
beam. 

2.4. Correction for CO emission 

Using the 19-channel bolometer, with its central frequency 
of 240 GHz and its effective bandwidth of Az/f, ; w 80 GHz, 
the 12 CO(2-l) line at 230 GHz might noticeably con- 
tribute to the total intensity in the bandpass. To get a 
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Fig. 1. Maps of M82 at A1.2 mm. The upper map shows 
the continuum emission after standard data reduction. 
The lower one shows the continuum emission, with the 
CO line emission subtracted. The contours are the same 
in both images starting at the 2er (10 mJy/beam) level. 
The beamsize of the IRAM 30-m telescope is shown in the 
lower left corner of each image. 



reliable information on the distribution of the cold dust, 
we applied the following correction to our continuum map: 
We used a 3' x 3' map of M 82 made with the same tele- 
scope in the 12 CO(2-l) line (obtained to provide the zero 
spacings for an investigation with the Plateau de Bure in- 



terferometer, Weifi et al.in prep.). This map was measured 
in the on-the-fly mode. The on-the-fly technique ensures 
a uniform sampling of the emission over the whole field 
with good signal-to-noise ratio. Since we knew from our 
continuum map that the effects of the error beam can be 
neglected at a frequency of 240 GHz, we did not apply any 
error beam correction to the CO map. 

Integrating over the effective bandwidth of the bolome- 
ter, the total continuum flux in the beam at the peak po- 
sition is given by 



S bo1 = A^ bo i • S^o = 2.6 • 10 



bol 
241 1 



-16 



Wm" 



where S||J = 330 mJy is the peak flux at 240 GHz. 

The integrated intensity of the 12 CO(2— 1) line at 
the same position jT m bdv = 625Kkms 1 was taken 
from our CO map. Expressing this in terms of fre- 
quency instead of velocity, the integrated intensity reads 
f T m b dv — 500 K MHz. 

To obtain the total flux in the 12 CO(2— 1) line we 
used the Rayleigh- Jeans Approximation and assumed 
that T m b ~ T cx : 



s co = Aj/co • n 



2kT n 



2kiv 2 



- mb 



With the beam solid angle £1 = 4.5 • 10 9 sr and 
the integrated 12 CO(2— 1) intensity from above we get 
S co =3.6-10- 17 Wm- 2 . Thus, S co /S bo1 « 0.14 at the 
position of the continuum peak. 

If we assume that this value holds for the whole galaxy, 
which of course is a simplification, we can obtain a CO- 
corrected continuum map by subtracting the appropri- 
ately scaled CO map from our continuum map. The un- 
corrected and the corrected map are shown in Fig. [l]. 

3. Results 

3.1. Overall distribution of the dust emission 

Fig. [2] shows the CO-corrected continuum map of M 82 
at A1.2mm, superimposed onto an optical image. At 
A1.2mm M82 has a more or less oval shape (ellipticity 
e f» 0.6), with the major axis parallel to the galactic plane 
(molecular gas). This corresponds to a position angle of 
approximately 70° . The maximum of the continuum emis- 
sion is located at a 50 = 9 h 51 m 41 s , A 50 = 69°54'55". 
Within the errors this position is identical to the peaks of 
the distributions at A800/iim and Al.lmm (Hughes etal. 
1990). Besides this regular shape there are some spur-like 
features extending above and below the galactic plane, 
which might be associated with the outflow seen at opti- 
cal and IR wavelengths. 

3.2. Major- and minor-axis profiles 

The contours in Fig. ^ indicate some extended emission 
also along the minor axis of the galaxy, but at first glance 
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Fig. 2. CO-corrected map of M82 
at A1.2 mm overlayed on a R-band 
image made with the 1.2-m tele- 
scope on Calar Alto. The contours 
start at the 2a (10 mJy/beam) level. 
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corner. 



the extent in this direction is not as large as suggested 
by the maps of Kuno & Matsuo (1997) and Alton etal. 
(1999). We have therefore computed the intensity as a 
function of distance from the galactic center along the 
minor (major) axis, averaging over some 30" in the major- 
axis (minor-axis) direction. These profiles are displayed in 
Fig. [| For comparison the 12 CO(2— 1) line intensity, the 
850 fim (Alton et al. 1999) and 20 cm continuum profiles 
are shown, too. 

The A 1.2 mm emission shows the most concentrated 
distribution, significantly smaller than that at 850 /im. Es- 
pecially the scale height along the minor axis is very small 
with respect to the other wavelengths and the molecular 
emission. 

Contrary to the expectation the molecular emission is 
the most extended component. Possible explanations for 
this phenomenon are discussed in Sect. 4.2. 



3. 3. Integrated flux density 

Integrating over the continuum map we obtain a total flux 
density of S240 = 1-91 ± 0.13 Jy. This is lower than the 
value measured by Kriigel et al. (1990a), but their flux 
density was much more uncertain, owing to the lower sen- 
sitivity of their 1-channel bolometer. 



Fig. |4| shows a radio-to-IR spectrum of M 82 with the 
individual values listed in Tab. [l]. Our new 240 GHz mea- 
surement is completely consistent with the slope of the 
spectrum at this frequency (indicated by the dotted line), 
whereas the value derived by Kuno & Matsuo (1997) at 
157 GHz lies much too high with respect to all adjacent 
data points. 

The CO-corrected flux density is significantly lower 
than the flux density in our original map. This demon- 
strates that a line correction for continuum measurements 
in the relevant bands is indispensible. 

The total flux can be turned into a gas mass via 



TvL = 



'240 



• D 2 



B(T d )- K ' 240 M d 



With a distance of D = 3.25 Mpc, and using the same 
parameters as Kriigel etal. (1990a), i.e. a dust tempera- 
ture of Td = 30 K (Chini et al. 1989) and a dust absorp- 
tion coefficient re 2 4o = K 240 ' m = 0-003 cm 2 g _1 (Kriigel 
etal. 1990b), we derive a gas mass of M g = 7.5 • 10 8 M Q . 
If we adopted the Draine & Lee value K240 = 0.005 cm 2 g _1 
(Draine & Lee 1984) the resulting gas mass would be some- 
what lower (M g = 4.5 ■ 1O 8 M ). 

In this calculation we assumed a gas-to-dust ratio of 
approximately M g /M<j w 100. The corresponding dust 
mass is M d = 7.5 • 10 6 M Q . 
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Fig. 4. Radio-to-IR spectrum of 
M82. Our new value is indicated by 
the filled circle, the CO-corrected 
value by the open circle, and the 
115 GHz point of Kuno & Matsuo 
1997 by the filled square. The slope 
of the spectrum in the submm wave- 
length range is indicated by the dot- 
ted line. The error bars of our data 
points are smaller than the used 
symbols. 



4. Discussion 

4--1- A dust halo? 

The A2mm map (Kuno & Matsuo 1997) shows contin- 
uum emission with a box-like shape, indicating a dust 
halo around M82. Although the emission at A 1.2 mm and 
A2 mm should be similar, we do not see any /*-evidence 
for a significant halo contribution to the cold dust. This 
difference could be due to the error beam of the Nobeyama 
45-m telescope. 

Effects of the error beam and the antenna pattern 
have been inv estigated and largely removed from our map 
(see Sect. 2J5), while this is not mentioned in case of the 
Nobeyama measurements. The total flux density in the 
Nobeyama map is 1.7 Jy, which is much more than ex- 
pected at this wavelength (see Fig. [|). This extra flux 
might be due to emission entering through sidelobes and 
the error beam. 

The 850 /zm map made by Alton etal. (1999) with the 
JCMT shows continuum emission at the 50 mJy level up 
to 40" above the galactic plane. This translates to 15 mJy 
at A1.2mm, corresponding to the second contour (3<r) in 
our map. Considering only the overall extent of the emis- 
sion, the JCMT map is consistent with ours. The different 
shapes might be due to problems with the baseline sub- 
traction, because the area mapped by Alton et al. is only 
slightly larger than the extent of the continuum emission. 



4-2. Comparison of dust and CO 

In non-active galaxies like NGC 4565 (Neininger et al. 
1996) and NGC 5907 (Dumke et al. 1997) the molecular 



line emission is more concentrated towards the central re- 
gion than the thermal dust emission. As can be seen from 
Fig. HI and [| the galaxy M82 exhibits the opposite be- 
haviour. The CO emission appears to be much more ex- 
tended than the continuum. Since the maps in Fig. ^| are 
made with the same telescope at the same frequency, the 
difference can not be due to different beam patterns and 
error beams. 

The dust is heated by the UV radiation field, which is 
concentrated towards the dense central region, because the 
UV photons can hardly escape without being absorbed. 
Only a small fraction of the UV photons makes its way out 
to kpc distances from the galactic plane, where it might 
contribute to the observed polarized radiation. 

The molecular gas can be excited by the radiation field, 
but also by low-energy cosmic rays and by soft X-ray 
photons (see e.g. Glassgold & Langer 1973), which have 
a much larger scale height than the UV photons (Shopbell 
& Bland-Hawthorn 1998). 

Fig. U shows that the synchrotron intensity at A20 cm 
drops faster with distance from the galaxy plane than the 
CO and even the dust emission. However, at low radio fre- 
quencies there may be a significant radio halo (low-energy 
cosmic rays have longer lifetimes), which is already indi- 
cated by the 327-MHz map of Reuter et al. (1992). Hence, 
there may be abundant cosmic rays for heating. Such cos- 
mic rays (E ^100 McV) remain, however, invisible in the 
radio window, since even in the strong (B » 10 ... 50 fiG) 
magnetic field of M82 they produce synchrotron radia- 
tion at frequencies below about 10 MHz. Soft X-rays are 
also seen far out of the plane of M 82 (e.g. Bregman et al. 
1995). These circumstances altogether provide the likely 
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Fig. 3. Profiles along the minor and major axis of M82. 
For each direction the A1.2 mm continuum emission (this 
paper), the integrated CO(2— 1) line intensity (this paper), 
the A850 ^m continuum emission (Alton et al. 1999) and 
the radio continuum emission at A20 cm (Reuter etal. 
1992) are shown. All data are smoothed to a beam of 
15" HPBW. 



heating sources for any molecular material transported out 
of the disk of M82 into the halo, rendering the CO "over- 
luminous". This is also indicated by the low conversion 
ratio of CO line intensity to molecular gas column den- 
sity, X co = N H2 /I C o derived by Smith ctal. (1991). The 
existence of a soft X-ray and a low-frequency radio halo 
around M 82 readily explains the larger extent of the CO 
emission as compared to the cold dust. 

5. Summary 

In this paper we presented a new high-sensitivity map 
of the continuum emission at A 1.2 mm in the starburst 



Table 1. Flux densities of M82 from radio through IR 
wavelengths. 
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Kellermann et al. 1969 
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Klein etal. 1988 
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Kellermann et al. 1969 
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Doherty et al. 1969 
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Klein etal. 1988 
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galaxy M82. To get a reliable map of the extended ther- 
mal dust emission we applied a CLEAN algorithm to our 
map (error beam correction) and subsequently subtracted 
the 12 CO(2-l) line emission. Besides this the contribution 
of the sidelobes of the IRAM 30-m telescope has been es- 
timated. 

The total mass of dust in the inner 3 kpc of the galaxy 
derived from the integrated A1.2 mm flux density is 5.0T0 6 
M Q , the inferred total mass of gas is 7.5 • 10 8 M Q . 
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Fig. 5. Comparison of the A1.2 mm continuum and the 
CO(2— 1) line emission. The upper image shows the CO- 
corrected A1.2 mm continuum map of M82. Contours 
start at 2cr (10 mJy/bcam). The lower image shows the 
integrated CO (2—1) line intensity. Contours start at 3a 
(20 K km s -1 /beam). The beam size of the TRAM 30-m 
telescope is given in the lower left corner of each image. 



Although the continuum emission is not confined to 
the galactic disk, our map rules out the existence of any 
pronounced halo of cold dust around M82. We find the 
CO(2— 1) line emission to be clearly more extended than 



the A1.2 mm dust emission. This unexpected observational 
result now waits for a conclusive theoretical explanation. 
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